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Abstract 
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The salient features of mixed-phase and ice clouds in a GCM cloud scheme are exam- 
ined using the ice formation parameterizations of Liu and Penner (LP) and Barahona 
and Nenes (BN). The performance of LP and BN ice nucleation parameterizations were 
5 assessed in the GEOS-5 AGCM using the McRAS-AC cloud microphysics framework 
in single column mode. Four dimensional assimilated data from the intensive observa- 
tion period of ARM TWP-ICE campaign was used to drive the fluxes and lateral forcing. 
Simulation experiments where established to test the impact of each parameterization 
in the resulting cloud fields. Three commonly used IN spectra were utilized in the BN 
10 parameterization to described the availability of IN for heterogeneous ice nucleation. 

The results show large similarities in the cirrus cloud regime between all the schemes 
tested, in which ice crystal concentrations were within a factor of 10 regardless of the 
parameterization used. In mixed-phase clouds there are some persistent differences 
in cloud particle number concentration and size, as well as in cloud fraction, ice water 
is mixing ratio, and ice water path. Contact freezing in the simulated mixed-phase clouds 
contributed to transfer liquid to ice efficiently, so that on average, the clouds were fully 
glaciated at T ~ 260 K, irrespective of the ice nucleation parameterization used. Com- 
parison of simulated ice water path to available satellite derived observations were also 
performed, finding that all the schemes tested with the BN parameterization predicted 
20 average values of IWP within ±15% of the observations. 

1 Introduction 

g 

__ C/)' 

The role of atmospheric aerosols in modulating the atmospheric radiative balance, by 
directly scattering solar radiation, or indirectly, modifying cloud optical and microphys- 
ical properties, has received considerable attention during the last couple of decades. 

25 Soluble and insoluble aerosol species provide nucleation sites for the atmospheric wa- 
ter vapor to form liquid droplets (Cloud Condensation Nuclei, CCN), and ice crystals 
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(Ice Nuclei, IN), respectively. The important interactions between aerosol particles and 
cloud optical and physical properties operate at temporal and spatial scales unresolved 
by Global Climate Models (GCMs); their inclusion in climate simulations therefore re- 
lies on parameterizations. The importance of these aerosol-cloud interactions, and 
5 their potential impact on climate, makes their inclusion in climate models through ac- 
curate and physically based schemes a high priority (Intergovernmental Panel on Cli- 
mate Change, 2007). 

Aerosol indirect effects (AIE) in warm clouds have been long studied and imple- 
mented in atmospheric models (e.g., Penner et al., 2006), but less has been accom- 
io plished for cold clouds. Modifications to the number density and sizes of ice crystals 
not only strongly affect the radiative properties of ice-bearing clouds, but also impacts 
the development of precipitation (e.g., Lohmann and Diehl, 2005; Lohmann, 2002). 

The complexities associated with cold and mixed-phase clouds (due in part to the con- 
current action of different freezing mechanisms, the high selectivity of the IN process, 

15 and the theoretical uncertainties associated with their description) have challenged the 
representation of such clouds in GCMs, most of which lack explicit ice microphysics 
(Lohmann and Feichter, 2005). As a result, even the sign of the radiative effects of 
aerosol-ice cloud interactions remains uncertain in climate simulations. 

(J) 

Important steps to improve the simple treatments of cold and mixed-phase cloud 
20 microphysics originally included in GCMs have been undertaken in recent years. For 
example, the partitioning of cloud condensate between ice and liquid water in mixed- 
phase clouds (235 K < T < 273 K) was typically represented by a temperature-only ap- 
proach (e.g., DelGenio et al., 1996; Rasch and Kristjansson, 1998). This approach has 
been progressively replaced by a less empirical and more physically-based represen- 
25 tation, in which the deposition growth of cloud ice at the expense of the liquid water, the 
Bergeron-Findeisen process (Pruppacher and Klett, 1997), is taken into account (e.g., 
Rotstayn et al., 2000). This prognostic approach for condensate partitioning which in- 
cludes explicit dependency of the deposition rate on microphysical variables such as 
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ice content q v and ice crystal concentration, A/ c , has been adopted by a variety of 
GCMs (Sud and Lee, 2007; Liu et al., 2007; Salzmann et al., 2010). 

Another advancement in GCM cloud schemes is the implementation of two-moment 
cloud microphyisics, which include prognostic equations for the mass as well as the the 
5 number concentration of different hydrometeor categories (e.g., Seifert and Beheng, 
2001, 2006; Morrison and Gettelman, 2008). This has permitted the prognostic com- 
putation of cloud particles sizes and ice deposition rate (e.g., Salzmann et al., 2010; 
Muhlbauer and Lohmann, 2009). 

Estimates of the AIE on ice-bearing clouds require an adequate description of the 
io aerosol-cloud coupling through the nucleation process. This is, the prognostic calcu- 
lation of hydrometeor sizes should be done in a manner consistent to aerosol load 
changes and aerosol characteristics. However, an efficient and comprehensive repre- 
sentation of the current understanding of the ice nucleation process in the framework 
of a GCM has proven difficult. Most ice nucleation parameterizations rely on simple 
15 functions to determine how many ice crystals will be heterogeneously nucleated at 
a given set of environmental conditions. These relations describing the availability of 
IN, termed IN spectra, exhibit different level of complexity, ranging from saturation- g 

dependent schemes (Meyers et al., 1992; Phillips et al., 2007) to IN spectra with 
aerosol-dependent parameters derived empirically (Phillips et al., 2008). Theory-based 
20 approaches have also lead to formulations of IN spectra with explicit dependence on 
aerosol number concentration, aerosol size distribution, and aerosol surface proper- 
ties (Khvorostyanov and Curry, 2009; Barahona and Nenes, 2009b; Barahona, 2012). 

Most GCM microphysical schemes that account explicitly for aerosol effects represent 
ice nucleation assuming that there is no variation in ice nucleation properties within 
25 an aerosol species. In reality, there is large variability in the ice nucleation properties 
of aerosol populations, which contributes to the large uncertainty in the predicted IN 
concentrations. 

Homogeneous freezing of solution droplets (i.e., without the presence of a solid 
aerosol phase) may occur only at temperatures below 235 K, the homogeneous 
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freezing threshold 7 hom (Pruppacher and Klett, 1997). For temperatures higher than 
7 hom , in which mixed-phase clouds typically exist, the presence of a solid phase is 
necessary for ice formation, and therefore only heterogeneous ice nucleation is ac- 
tive. Below 7 hom , where ice-only clouds form, the supersaturation with respect to ice is 
5 the result of the competition between the rate of cooling of the cloud parcel and the 
condensation on the nucleated ice crystals. Therefore, it varies dynamically given the 
amount of IN present and the dynamical forcing available. Furthermore, since homo- 
geneous and heterogeneous ice nucleation may occur simultaneously, the competition 
from both mechanisms and their impact on supersaturation further complicate the cal- 
10 culations. For this reason Lagrangian simulations have been used to develop solutions 
to the variable supersaturation problem (e.g., Lin et al., 2002), and fits to these nu- 
merical solutions have been used to develop ice nucleation parameterizations. A few 
such parameterizations have been developed (Karcher and Lohmann, 2002, 2003; Liu 
and Penner, 2005), and have been implemented in GCM models (Hoose et al., 2010). 
15 Analytical solutions to this problem have been developed in which any IN spectra can 
be used (Barahona and Nenes, 2009b). 

For the case of mixed-phase clouds, liquid water, water vapor, and ice are simultane- 
ously present, and can exhibit complex dynamics (e.g., Korolev, 2007). For the coarse 
resolution of GCM cloud schemes, the simplifying assumption that the water vapor is 
20 saturated with respect to liquid water is sometimes made. The supersaturation with 
respect to ice, S-„ is therefore constrained by thermodynamic equilibrium rather than 
by the competition of cooling and condensation. With this assumption, it is sufficient to 
know the availability of IN (given by an IN spectrum) at S, to compute the nucleation 
rate of ice crystals. 

25 A number of studies have focused on the implementation and evaluation of new 
microphysical schemes in GCM simulations, including prognostic calculation of the 
ice fraction in mixed phase clouds, and using more physically-based ice nucleation 
schemes (Storelvmo et al., 2008; Sud and Lee, 2007; Liu et al., 2007; Salzmann 
et al., 2010). Curry and Khvorostyanov (2012) performed a comparison of some 
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heterogeneous nucleations parameterizations in a single-column model for long lived 
mixed-phase arctic clouds. However, none of these studies have performed sensitivity 
analysis of the simulated mixed-phase and cirrus clouds fields to different IN spectra. 

In this study, we use the parameterization of Barahona and Nenes (2009b), BN here- 
5 after, in which the ice nucleation problem is treated in a general framework that admits 
the use of any IN spectra, empirical or theoretical. Barahona et al. (2010) used the BN 
parameterization to compare common formulations of the IN spectrum in a chemical 
transport model, finding that the 2 to 3 orders of magnitude variation in the IN con- 
centrations among different schemes would lead to up to a factor of 20 variation in N c 
io in cirrus clouds. The sensitivity can be even larger in mixed-phase clouds where only 
heterogeneous ice nucleation is active, and competition for water vapor does not buffer 
the response of crystal number to IN concentration changes. 

Testing the impact of IN spectra in a comprehensive cloud microphysical framework 
would provide valuable information on how the uncertainties associated with ice nucle- 
15 ation are reflected on the cloud field variables when coupled to other cloud processes. 
In this study, we report the implementation of the BN ice nucleation scheme into the 
Microphysics of Clouds with Relaxed Arakawa-Schubert and Aerosol-Cloud interaction 
(McRAS-AC) (Sud and Lee, 2007) driven by the Goddard Earth Observing System 
Model, version 5 (GEOS-5). The flexibility provided by the BN ice nucleation parame- 
20 terization is ideal for testing the sensitivity of the simulated cloud properties to the rep- 
resentation of IN spectra in the McRAS-AC framework. To isolate the response from the 
underlying physical parameterization, all the simulations were performed in the Single 
Column Model version of GEOS-5. This is a common test of GCM microphysics since 
the SCM configuration contains the same physical parameterizations as the host GCM 
25 model, with the advantage of a much smaller computational burden, and the laterally 
constrained input flux fields allows better delineation of the role of microphysical pro- 
cesses of cloud formation and aerosol effects. The simulations were forced with data 
collected during the Tropical Warm Pool International Cloud Experiment (TWP-ICE) 
intensive observation period (IOP) of the ARM program (May et al., 2008). 
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2 Model description and simulation set-up 

A detailed description of the McRAS-AC microphysics can be found elsewhere (Sud 
and Walker, 1999; Sud and Lee, 2007; Bhattacharjee et al., 2010). Here we will pri- 
marily focus on describing the treatment of cold and mixed-phase clouds microphysics 
5 in McRAS-AC. 

2.1 Ice nucleation in McRAS-AC 

McRAS-AC has the option to invoke the Liu and Penner (2005), (LP), or the Bara- 
hona and Nenes (2008, 2009a, b) parameterizations to describe the ice nucleation 
process. This possibility was used to assess and compare the performance of the 
io two schemes. The LP parameterization was originally designed to describe the nucle- 
ation process at temperatures typical of cirrus cloud formation, i.e., for temperature 
less than the homogeneous freezing threshold (T hom = 235 K). It is based on numerical 
correlations derived from statistical fits to a large number of Lagrangian parcel model 
simulations, in which homogeneous and heterogeneous freezing mechanisms were 
15 explicitly accounted for. The homogeneous freezing of deliquesced sulfate aerosol was 
approached using an effective freezing temperature. Immersion freezing on soot parti- 
cles was included in the parcel model simulations using a classical nucleation theory 
description, in which a fixed aerosol size distribution and freezing characteristics were 
assumed. Deposition freezing is calculated using the Meyers et al. (1992) formulation. 
20 In this way, the LP parameterization takes into consideration the impact of updraft ve- 
locity, w , and aerosol load N a on the number concentration of nucleated ice crystals, 
N c nuc . In the cirrus regime, it is capable of calculating N c nuc as the result of the com- 
petition of both freezing mechanisms. For temperatures above 7" hom , homogeneous 
freezing is inactive, and A/ cnuc is given solely by heterogeneous nucleation. However, 
25 because the scheme obtains the number of nucleated ice crystals from curve-fitted 
functions of temperature and vertical velocity, these specific equations may not hold as 
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well when extrapolated beyond the curve-fit data domain, nor be applied to aerosols 
that does not follow the prescribed freezing properties used on the simulations. 

The implementation of LP in McRAS-AC for mixed-phase clouds (7" hom < T < 273 K) 
follows closely that of Liu et al. (2007), and was described and tested in a SCM frame- 
work (Bhattacharjee et al., 2010). In this regime, A/ cnuc is calculated by adding the 
contributions from the the numerical correlations described above and the contribution 
from deposition freezing as given by a modified version of the Meyers et al. (1992) 
formula, 

o 

N\d(S\) = f (z)N 0 exp(a + b(S t - 1 )) (1) | 

(J) 

where A/ id is the number concentration of ice crystals due to deposition nucleation 
in m -3 , N 0 = 10 _3 m“ 3 , a = -0.639, and b = 0.1296, and f(z) is an empirical height 
correction factor, given by f(z ) = io (z ° -z)/, ' 5z , with z 0 = 1 km, Sz = 6.7km, and f{z) e 
[0.12,1.0]. This decay factor was derived from observations by Minikin et al. (2003) 
during the INCA (Interhemispheric Differences in Cirrus Properties from Anthropogenic 
Emissions) campaign, to augment the formula by Meyers et al. (1 992) that was derived 
from ground-level observations. 

In the present work we implemented and tested the BN parameterization in McRAS- 
AC. BN is based on an analytical solution of the governing equations of a cooling 
air parcel in which deliquesced aerosol and heterogeneous IN are allowed to freeze 
and grow by water vapor deposition (Barahona and Nenes, 2008, 2009a, b). Accord- 
ingly, BN circumvents the need for curve-fitted equations, and holds for a wide range 
of configurations encountered in the physical system. The availability of IN in the BN 
parameterization can be described, in principle, with any heterogeneous nucleation pa- 
rameterizations. Here we use the correlations of Meyers et al. (1992) (MY92), Phillips 
et al. (2008) (PDA08), and the semi-empirical spectra derived from classical nucleation 
theory of Barahona and Nenes (2009b), (CNT BN). MY92 is a widely used, empirical 
IN spectrum which depends only on S r PDA08 is also empirical however it separately 
considers the contribution of organics, dust, and black carbon to the IN population. 

14935 


ACPD 

12, 14927-14957, 2012 


Sensitivity of 
mixed-phase clouds 
to IN spectrum 

R. Morales Betancourt 
et al. 


Title Page 




Full Screen / Esc 


Printer-friendly Version 


Interactive Discussion 






CNT BN is based on an approximation of classical nucleation theory to scale observed 
IN concentration as a function of supersaturation. For temperatures T < 7" hom , BN calcu- 
lates the competing effects of homogeneous nucleation on deliquesced aerosol and the 
heterogeneous freezing for the availability of water vapor in a forming cirrus cloud. The 

5 maximum supersaturation with respect to ice attained in the ascending parcel, S- t max , is 
calculated by balancing the depletion effect from deposition growth of ice crystals and 
the availability of water vapor from cooling. In this way, S i>max is given by the dynamics 
of cooling and ice nucleation. BN then uses the maximum saturation to calculate N c nuc . 

The application of BN in the mixed-phase cloud regime differs slightly from that of LP. 

10 In the absence of any liquid water, the maximum supersaturation in the parcel would 
be dictated dynamically by expansion cooling and by the IN concentration. However, in 
McRAS-AC, any initial condensate is considered to be liquid (Rotstayn, 1997) and its 
then partitioned following Rotstayn et al. (2000). Therefore, in practice, ice nucleation 
above 7 hom is assumed to occur in an environment saturated with respect to water. 

15 Under this circumstance, S, max is fixed by the assumption of water saturation, equal 
to S i max = e sw (7")/e si (7"), i.e. , the ratio of the saturation vapor pressure over water and 
over ice, and it is therefore independent of the dynamic forcing, w, or aerosol loading. 

The number concentration of nucleated ice crystals, N c nuc , is then calculated by direct 
application of the IN spectra at the given S t max . 

“O 

20 2.2 McRAS-AC cold and mixed-phase cloud microphysics 

— i 

The cloud microphysics in McRAS-AC include balance equations for the mixing ratios 
of liquid water, q x , and cloud ice q r The precipitation microphysics are described by 
Sud and Lee (2007), which recast Seifert and Beheng (2006) to apply it to the thicker 
clouds of a coarse resolution GCM. The activation of aerosol to cloud droplets follows 

C D 

25 the parameterization of Fountoukis and Nenes (2005). Aerosol mass concentrations 
are taken from the Goddard Chemistry Aerosol Radiation and Transport (GOCART), 
and log-normal size distribution for each species are prescribed. 
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The partitioning of cloud condensate between ice and liquid in mixed-phase clouds 
is prognostic, and takes into account the Bergeron-Findeisen (BF) process, by which 
cloud droplets evaporate and the resulting water vapor deposits, to ice crystals. The 
process is represented following Rotstayn et al. (2000), which explicitly accounts for 
the dependence of the ice deposition rate on crystal number concentration, N c . 

Ice crystal number concentration N c is determined in McRAS-AC by the processes of 
ice nucleation, contact freezing, and by melting of cloud ice. The ice nucleation term is 
calculated with the LP and the BN parameterizations as explained in Sect. 2.1 . Contact 
freezing of supercooled cloud droplets through Brownian coagulation with insoluble IN 
(mineral dust) is included as given by Young (1974). 

Aerosol input for ice nucleation is also based on GOCART aerosol climatology. A sin- 
gle mode log-normal size distribution was assumed for black carbon, with geometrical 
mean diameter, d g = 0.04pm, and a geometric standard deviation o g = 2.3 (Jensen 
and Toon, 1994). Similarly, sulfate aerosol size distribution is assumed log-normal with 
d g = 0.14pm and o g = 1.5 (Pueschel et al., 1992). The density of black carbon was 

assumed equal to 1 gem -3 while for sulfates, we assumed the density of sulfuric acid 
(1.84 gem -3 ). A probability distribution function of cloud scale vertical velocity, w, was 
used to represent the local variations of velocity at scales relevant for nucleation. The 
distribution was assumed to be a normal distribution with a fixed standard deviation 
of 0.25 ms -1 , consistent with observations in the INCA campaign (Karcher and Strom, 
2003). 

2.3 Forcing data 

g 

The SCM configuration consists of an isolated column of a global circulation model, 
and is therefore, a 1 -dimensional time-dependent atmospheric model. The lateral forc- 
ing fields to the 72 pressure levels in the atmospheric columns of GEOS-5 are pre- 
scribed from assimilated 4-D observational data. For the purpose of this study, we 
used the forcing from the TWP-ICE intensive observation period (IOP), derived by the 

- \ 
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Atmospheric Radiation Measurement (ARM) program. It includes data from 1 7 January 
to 12 February 2006. This data set has been previously utilized in forcing SCM sim- 
ulations with the intent of testing ice microphysics for GCMs (Wang et al., 2009a), as 
well as for comparing simulations produced with bulk microphysical schemes of varying 
complexity in a cloud resolving model with observational data (Wang et al., 2009b; Lee 
and Donner, 2011). The TWP-ICE data is ideally suited for testing the representation 
of cold and mixed-phase clouds in models and is an often used test case that allows 
comparison with other existing studies (e.g., Varble et al., 2011; Fridlind et al., 2012). 

It includes periods dominated by deep convective clouds and by persisting layers of 
cirrus clouds. 

(J) 

O' 

15 

3 Simulated clouds fields 

"0 

CD 

—i 

Table 1 summarizes the simulations considered in this study. The objective of the sim- 
ulation experiments is to evaluate the sensitivity of the cloud fields to the treatment 
of ice nucleation. A “control” simulation was performed with the LP ice nucleation pa- 
rameterization as described above, which has been used in the McRAS-AC framework 
before (Bhattacharjee et al., 201 0). Other simulation experiments were carried out with 
the BN parameterization, utilizing three different IN spectra. Two additional simulations 
were considered, in which the contribution to N c from contact freezing was neglected 
(LP-NoFrzc and BN-PDA08-NoFrzc). All the simulations share the same lateral forcing 
fields, surface fluxes, and aerosol input, and they only differ on the treatment of ice 
formation. 

The time-height distributions of the total cloud fraction, CF, exhibits the basic fea- 
tures observed during the TWP-ICE campaign (Fig. 1). In the first period of the in- 
tensive observation period (IOP), prior to 25 January 2006, the region was influenced 
by an active monsoon period characterized by considerable convective activity. From 
26 January to 2 February, the monsoon was suppressed, and little convective activity 
was observed, but high clouds persisted through the period. In the final part of the IOP 
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(3-13 February) the region was increasingly impacted by continental storms, reflected 
in a renewed increase in the convective activity. 

The simulated cloud fields show some differences in the CF, particularly the simula- 
tion with the BN-PDA08, which shows higher frequency of high CF cells. Common to 
5 all the simulations with the BN scheme is an increase in the CF for the mixed-phase 
regime as compared with the LP-CTRL simulation, particularly in the convectively ac- 
tive periods, as shown for two of the simulations in Fig. 3. The resulting simulated q t 
fields are shown in Fig. 2. Ice mixing ratios generally reach a maxima in the layer ex- 
tending from the 0°C to the -38 °C levels. The overall ice mixing ratios encountered 
io in the LP-CTRL simulation are generally higher than for the BN cases, the difference 
being more pronounced for the mixed-phase regime. 

The temperature dependence of A/ cnuc and the total ice crystal concentration N c 
was calculated from the model output for each one of the simulations as a function 
of temperature (Fig. 4). The impact of the nucleation scheme in the partitioning of 
is condensate was investigated through the ice fraction, f c , defined as 


Figure 5 shows the temperature dependence of the condensate partitioning, for the 
BN-PDA08, LP-CTRL, BN-PDA08-NoFrzc, and LP-NoFrzc. Attention was given to vari- 
ables affecting the radiative properties of the ice clouds. The size of ice particles would 
20 be among the most directly affected variables with changes in crystal concentrations. 
The behavior of effective radius for ice particles as a function of temperature is shown 
in Fig. 6 for two of the simulations. Figure 7 shows a time series of IWP from differ- 
ent simulation experiments, together with IWP derived from satellite retrievals using 
the Visible Infrared Shortwave-Infrared Split-Window Techinique (VISST), described in 
25 Fridlind et al. (2012). 
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4 Discussion of the results 

Since the lateral forcing and surface fluxes were prescribed identically in all simulation 
experiments, any differences in the cloud fields can be attributed to the interaction of 
the ice nucleation scheme with the BF process and the cloud microphysical response 
5 that follows. Some such differences are encountered between the fields produced with 
LP and BN parameterization, respectively. /V cnuc calculated with BN is systematically 
lower for the mixed-phase cloud regime irrespective of the heterogeneous nucleation 
scheme used, however, the difference is greater between LP-CTRL and BN-PDA08, 
for which the maximum difference in the predicted A/ cnuc can be considerable (Fig. 4). 
io The low concentration of IN predicted by the PDA08 spectrum, typically two orders of 
magnitude lower than produced by the other spectra, explains part of this difference. 
Flowever, the systematic discrepancy between LP and BN in the mixed-phase regime is 
likely due to the different implementation of the two nucleation schemes. As described 
in Sect. 2.1, the LP scheme adds the contributions from immersion freezing (given by 
is the numerical correlations of Liu and Penner, 2005, and from deposition, as given by 
Eq. 1). In the BN schemes the availability of IN in the mixed-phase regime is dictated 
by the IN spectrum alone, which consider deposition and condensation freezing. 

The large differences in predicted N c nuc are also noticeable in the resulting N c fields, 
but the magnitude of the difference is significantly lower. In the range of temperatures 
20 where contact freezing is active (270.1 5K > T > 235 K) this mechanism was found to 
contribute, on average, between 10 _ 4 cm “ 3 and 10 _ 3 cm “ 3 to the ice crystal concen- 
tration, thereby effectively providing a lower bound for A/ c (Fig. 4). This contribution is 
significant only for IN spectrum predicting very low /V cnuc (such as PDA08), or for the 
temperatures above T ~260K, in which the other IN spectrum (MY92 and BN-CNT) 
25 predicts very small N c nuc . 

It is expected that the differences in N c would significantly impact other cloud micro- 
physical variables, particularly through the modification of the rate of the BF process. 
Lower ice crystal concentrations should result in lower rates of conversion of liquid 
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water to ice because the surface area for vapor-ice mass transfer is low (Rotstayn 
et al., 2000). Such behavior, in which low aerosol concentrations are associated with 
low f c , has been observed in satellite retrievals (Choi et al., 2010). However, the ice 
fraction exhibits little to no change across simulations even for the cases where N c 
differ by a factor of 100 (Fig. 5a, b). This diminished sensitivity of f c to ice crystal con- 
centration seems to be caused by the action of the contact freezing mechanism. To 
verify this, two simulations in which this mechanism was neglected were performed 
LP-NoFrzc and BN-PDA08-NoFrzc, (Fig. 5c, d). LP-NoFrzc shows that the transition 
from pure liquid to pure ice cloud occurs over a larger temperature interval as com- 
pared to simulations in which contact freezing is allowed to occur. However, because 
LP predicts relatively large crystal concentrations in the entire range of supercooling 
temperatures, the BF process is always fast, resulting in a rather similar dependence 
of f c on temperature. This is not the case for the simulations with BN-PDA08, in which 
the low N c severely limits the rate of conversion of liquid water to ice by water vapor 
deposition, which is evidenced when contact freezing is turned off (Fig. 5d). The nar- 
row temperature range associated with the transition from f c = 0 at 273 K to f c = 1 at 
260 K when contact freezing was included, is consistent with other studies with the 
same partitioning scheme, as well as with available cloud observations of f c (Rotstayn 
et al., 2000; Liu et al., 2007). The transition to a fully glaciated state which sometimes 
produces excessive ice cloud at low temperatures could arise from the assumption that 
the BF mechanism dominates this regime (Korolev, 2007). 

The cloud amount in the mixed-phase regime was affected by the crystal ice nucle- 
ation parameterization used in the simulations. As shown in Fig. 3, there is an increase 
in the frequency of occurrence of cloudy cells with CF > 0.5 when the BN-PDA08 pa- 
rameterization is used instead of LP. This is true for the three IN spectra utilized in this 
study, with CF being 49% larger for BN-PDA08, and ~ 25% for MY92 and BN-CNT, as 
compared to simulations with LP. 

In the cirrus cloud regime, the difference in A/ cnuc between LP and BN is less 
pronounced than in the mixed-phase regime. For this temperature range, crystal 
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concentrations calculated with LP and BN are within one order of magnitude irrespec- 
tive of the IN spectrum used, which is consistent with the variability reported in previous 
studies (Barahona et al., 2010). A/ cnuc for LP-CTRL and BN-PDA08 are in close agree- 
ment, specially for temperatures above 200 K, however, the predicted mechanism of 
5 freezing is different for both parameterizations. Due to the very low IN number pre- 
dicted with PDA08, the contribution of heterogeneous freezing to A/ c nuc in BN-PDA08 
is negligible, and the process is dominated by homogeneous freezing. The opposite 
behavior is observed when the LP parameterization is used, in which homogeneous 
freezing only contributes significantly to the N c nuc at extremely low temperatures. When 
io BN-CNT or MY92 are used instead, the lower N c nuc is the result of the depletion of wa- 
ter vapor from the more numerous IN, and homogeneous freezing is triggered only at 
temperatures between 200 K and 220 K (Fig. 4). 

Finally, even though the impact of N c on the simulated condensate partitioning is 
small, the different ice crystal concentration predicted with the parameterizations con- 
15 sidered in this study considerably impact the cloud radiative properties. For instance, 
Figs. 2 and 3 shows the ice mixing ratio for different simulation scenarios. The differ- 
ences observed translate also into ice water path differences, as well as of the hy- 
drometeor sizes. Figure 6 shows the temperature dependence of the median values 
of the calculated ice effective radius for BN-PDA08 and LP-CTRL. The inset shows 
20 a histogram of the frequency distribution of the effective radius for BN-PDA08 and LP 
only for the rage of mixed-phase temperatures. Due to the much lower N c predicted 
by PDA08, the effective radius is shifted from a median of 45 pm for BN-PDA08, to 
a smaller size with a median of 32 pm in the LP-CTRL simulation. 

The changes induced in the cloud microphysics by the different IN spectra conse- 
25 quently modify the overall column integrated properties of the cloud fields. Figure 7 
illustrate the tendencies in the simulated IWP for the different parameterizations. It is 
clear that IWP for LP-CTRL is higher than for any simulation with BN, with differences 
being larger in the periods of convective activity. In fact, the average IWP in the ac- 
tive monsoon period for the LP simulation was found to be 0.30 kg m -2 , while it was 
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of 0.22 kg m -2 for BN-PDA08, and 0.28 kg m -2 for BN-CNT and BN-MY92. In the sup- 
pressed period, IWP averaged ~ 0.04 kg m -2 in all the simulation experiments. This 
simulation results compare qualitatively well to the available data of IWP as retrieved 
from VISST. However, the lower bound in the VISST observations tends to be much 
lower than simulated IWP, while the peaks during the convective events often exhibit 
higher values than simulated fields. Ice water path from VISST during the active mon- 
soon period averages 0.25 kg rrf 2 and 0.04 kg rrf 2 for the suppressed period. 

D 

</)’ 

5 Summary and conclusions 

C n 

o' 

The ice nucleation parameterization of Barahona and Nenes (2009b) was implemented 
in the GEOS-5 McRAS-AC cloud scheme and tested in single column mode forced 
with TWP-ICE campaign data. Three different heterogeneous ice nucleation spec- 
tra (PDA08, BN-CNT, and MY92) were used in simulations experiments with the BN 
parameterization framework. The IN concentration predicted by the spectra used in 
this study varied greatly, with PDA08 predicting very low IN concentrations of around 
~10 _4 cm“ 3 , while MY92 and BN-CNT predicted IN concentrations similar to each 
other, but generally ~ 100 times larger than PDA08 at any given temperature. These 
simulation experiments were compared to a control simulation using the LP parame- 
terization, which was found to predict the highest ice crystal concentrations across the 
simulations. 

It was shown that the different schemes used in this study often predicted IN con- 
centrations differing by up to three orders of magnitude. Despite these important differ- 
ences in IN availability, ice crystal number concentration for cirrus cloud temperatures 
predicted in all the simulations were found to agree within a factor of 10. However, 
the mechanism by which these ice crystal are produced is considerably different; N c 
computed with LP was dominated by heterogeneous freezing, while simulations with 
BN transitioned from heterogeneous to homogeneous dominated freezing at higher 
temperatures. 
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In the regime of mixed-phase clouds, the variations in N c among simulations with the 
different nucleation schemes was considerably larger than for the ice-only clouds, with 
the largest variations being within a factor of ~ 100 in some cases. This larger vari- 
ability is not surprising, since in the absence of homogeneous freezing, the nucleation 
5 schemes strongly depend on the IN nucleation spectra. However, the contribution to 

o o 

N c from contact freezing of cloud droplets with dust particles of ~ 10 cm provided 
a lower bound on N c , and was effectively the largest contributor to crystal concentration 
when the PDA08 scheme was used. This contribution to N r also acted to counteract 
the very large variations in predicted IN concentrations. 
io Similarly, it was also found that the action of contact freezing efficiently transforming 
cloud water into cloud ice buffered the impact of the large variations of N c seen across 
the different simulation experiments on the partitioning of cloud condensate. Ice mix- 
ing ratios, however, where strongly affected by the ice nucleation scheme. Accordingly, 
cloud microphysical variables relevant to radiative properties, such as the effective ra- 
15 dius of ice crystals and the ice water path, were impacted by the wide range of N c 
predicted. It was observed that nucleation schemes that predict lower N c lead to lower 
in-cloud ice mixing ratios and ice water path, and considerably larger crystal sizes. 

This study highlights the need for detailed cloud microphysical observations to con- 
strain the large uncertainties associated with the ice nucleation process which limit the 
20 ability of GCM models to make accurate estimates of the contribution of cold clouds to 
the overall aerosol indirect effects. Continued development and refinement of ice nu- 
cleation schemes capable of accounting correctly for different freezing mechanisms is 
needed; using the approaches used here will help accomplish this. 

g 
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Table 1. Simulations reported in this study. The 
identical, and they only differ on the ice nucleation 


simulation 

scheme. 


set-up for all the simulations is 


Simulation ID 

Contact Freezing Ice-only Clouds (cirrus) 

Mixed-Phase Clouds 

LP-CTRL 

Yes 

Liu and Penner (2005) 

Meyers et al. (1992) 

LP-NoFrzc 

No 

Liu and Penner (2005) 

Meyers et al. (1992) 

BN-PDA08 

Yes 

Barahona and Nenes (2009a) 

Phillips et al. (2008) 

BN-PDA08-NoFrzc 

No 

Barahona and Nenes (2009a) 

Phillips et al. (2008) 

BN-CNT 

Yes 

Barahona and Nenes (2009a) 

Barahona and Nenes (2009a) 

BN-MY92 

Yes 

Barahona and Nenes (2009a) 

Meyers et al. (1992) 
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Fig. 1. Time-Height distribution of the simulated cloud fraction for (a) control simulation with 
the LP ice nucleation parameterization, (b) Simulation with BN and the PDA08 ice nucleation 
spectra, (c) Simulation BN and the CNT ice nucleation spectra, (d) Simulation BN09 with the 
MY92 ice nucleation spectra. The gray contours correspond to the 0°C and -38°C isolines, 
indicating the region where mixed-phase clouds may occur. 
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Fig. 3. (a) Histogram of cloud fraction for the LP-CTRL and BN-PDA08 simulations. The fre- 
quencies are calculated for the active monsoon period and for cells with temperatures in the 
range 235 K> T >273K. (b) Vertical profile of in-cloud ice mixing ratios averaged over the 
monsoon active period. 
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Fig. 4. (a) Average number of nucleated ice crystals, A/ cnuc as a function of temperature for 
the simulations considered in this study. The inset is the fraction of crystals nucleated het- 
erogeneously, A/ het //V c nuc . (b) Average number concentration of ice crystals N c as a function of 
temperature. The vertical dashed line marks the homogeneous freezing temperature threshold. 
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Fig. 5. Frequency distribution of the ice fraction f c as a function of temperature. The dark gray 
lines represent the quartiles of the distribution of f c corresponding to a temperature interval of 
1 K. (a) LP-CTRL simulation, (b) LP-NoFrzc simulation, (c) BN-PDA08 simulation, and (d) the 
BN-PDA08-NoFrzc. 
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Fig. 6. Median values for the ice crystals effective radius for BN-PDA08 and LP-CTRL. The 
inset is a histogram of the frequency distribution of the effective radius of ice crystals for the 
simulated clouds in the mixed-phase temperature regime (235K<7 <273K). The bins are 
uniformly separated in logarithmic scale. 
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Fig. 7. Ice Water Path (IWP) in kgm 2 from VISST data, and simulated IWP for different IN 
spectra (a) LP-CTRL and BN-PDA08, (b) LP-CTRL and BN-CNT, (c) LP-CTRL and BN-MY92. 
The dashed vertical lines denote the initiation and end of the suppressed monsoon period. 
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